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ABSTRACT 

We report the detection of ionized intracluster gas in the globular cluster 47 Tucanae. Pulsars in this cluster with 
a negative period derivative, which must lie in the distant half of the cluster, have significantly higher measured 
integrated electron column densities than the pulsars with a positive period derivative. We derive the plasma 
density within the central few pc of the cluster using two different methods which yield consistent values. Our 
best estimate of Wg = 0.067 ±0.015 cm"-^ is about 100 times the free electron density of the ISM in the vicinity of 
47 Tucanae, and the ionized gas is probably the dominant component of the intracluster medium. 

Subject headings: globular clusters: individual (47 Tucanae) — pulsars: general 



1. INTRODUCTION 

Every ~ 10^ yr, globular clusters pass through the plane of 
the Galaxy, during which any intracluster gas is expected to 
be stripped from the systems (Spergel 1991). However, winds 
from evolved stars continuously fill the cluster with gas be- 
tween passages. Heretofore, numerous searches for the ex- 
pected neutra l and ionized material have been unsu ccessful (see 
refer ences in [Spergel 199 1|; Cesser & Shawl 1977|; Smith et al 



1990; [Knapp et al. 1996| ; [Penny, Evans, & Odenkirchen 1997 ; 
Hopwood et al. 19991) or inconclusive (K rockenberger & Grind 



lay 1995; priglia et al. 1997b[ ). The most convincing evidence 
presented so far was a 3 cr detection of two CO lines in the direc- 
tion of 47 Tucanae (henceforth 47 Tuc). These were interpreted 
as resulting from a bow shock formed as the cluster moves in 
the Galactic halo (Origlia et al. 1997a), indicating the presence 
of an intracluster medium. There is also an unconfirmed 3.5 cr 
SCUBA detectio n of dust in NGC 6356 at A = 850 nm (H op- 
wood et al. 1998), and a tentative detection of diffuse X-ray 
emission by the intracluster medium of NGC 6779 (H opwood 
et al. 2000). 

In this Letter we report the detection of ionized gas in 47 Tuc 
from new measurements of the radio disp ersion of 15 of the 20 
millisecond pulsars known in the cluster ([Camilo et al. 200C ). 



2. OBSERVATIONS 

Observations of 47 Tuc using the Parkes radio telescope at 
a frequency of 1.4 GHz have been used to measure the pre- 
cise positions, periods P, and apparent rat es of change of pe - 
riod fobs, for 15 of its 20 known pulsars ( |^reire et al. 2001 ). 
We have continued timing observations of the pulsars in this 
cluster, and since 1999 August have used a new system pro- 
viding a three-fold increase in time resolution. This consists of 
the central beam of the Parkes telescope multibeam receiver, as 

'http://pulsar.princeton.edu/tempo 



before, but uses a 2 x 512 x 0.5-MHz filter bank centered on 
1390 MHz, and a sampling interval of 80 /is. With these data 
we acquire very high quality pulse profiles, from which we ob- 
tain the precise dispersion measure for each pulsar (DM, the 
integrated electron column density along the line of sight from 
the Earth), from the relative delays in the arrival times of pulses 
in four contiguous 64-MHz subbands. 

We do this with the TEMPO timing software' as follows. First 
we determine timing solutions for all pulsars including spin, as- 
trometric, and binary parameters where relevant, as detailed in 
Freire et al. (2001), but using only data collected at 1.4GHz 
(medium-resolution data during 1997 August-1999 August and 
high-resolution data during 1999 August-2001 February). In 
virtually all cases the parameters obtained are consistent at the 
3 (T level with those reported by Freire et al. (2001). We then 
hold the pulsar ephemerides fixed at these values, and fit for 
DM, using only high-resolution times-of-arrival measured in 
four frequency subbands. The resulting DMs are listed in Ta- 
ble nl along with the (P/P)ohs, and angular offset from the cen- 
ter of the cluster for each pulsar, 9±, as determined by Freire et 
al. (2001). 

As Figure |l] shows, all 15 pulsars are located near the cen- 
ter of the cluster Having periods of a few milliseconds, they 
are expected to have the small positive intrinsic period deriva- 
tives Pint which are typical of millisecond pulsars. However, the 
gravitational field of the cluster causes the pulsars to accelerate 
toward its center The line-of-sight component a of this acceler- 
ation results in a contribution to the observed period derivative 
of each pulsar given by P/P = a/c, where c is the speed of 
light. In the central regions of 47 Tuc, these values are typi- 
cally gr eater than the {P/P)-„, for millisecond pulsars (F reire et 
al. 2001), so that the pulsars act as tracers of the cluster gravi- 
tational field. 
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Fig. 1. — The positions of 15 millisecond pulsars in 47 Tuc superimposed 
on an optical image of the cluster (obtained by D. Malin, Anglo-Australian 
Observatory). All pulsars are loc ated within I '7 Ty 2 pc, delineated by the 
circle) of the center of the cluster ( Freire et al. 2001 ). North is to the top and 
east is to the left. 



Conversely, from a knowledge of the gravitational field of 
the globular cluster, one can use an observed period derivative 
to constrain the position of a pulsar in the cluster along the line 
of sight. In particular, pulsars located on the far side of the clus- 
ter will experience an acceleration directed towards the Earth, 
resulting in a negative period derivative which often exceeds 
the intrinsic value, leading to a negative value of iP/P)nbs- We 
observe this for nine of the 15 pulsars (Table p. Figure ^jshows 
the measured DMs plotted against (P/P)obs- Strikingly, most of 
the nine pulsars with negative (P/P)obs, which must lie in the 
distant half of the cluster, have significantly higher DMs than 
those with positive (P/P)obs, most of which are likely to reside 
on the near side. 

We now show that this configuration is unlikely to arise by 
chance, by computing the chance probability of observing (as 
we do) the 7 pulsars with the highest DMs all having a negative 
P (we neglect in this calculation 47 Tuc L and T, whose DMs 
and uncertainties make their contribution to this argument un- 
certain). The number of sets of 7 pulsars among the 13 with a 
timing solution is 13!/(7! x6!)= 1716. The number of sets of 7 
pulsars among the 8 with a negative P is 8. Therefore the prob- 
ability that the observed arrangement is due to chance is small, 
8/1716 « 0.005. It is conceivable that the observed variations 
in DM are simply a result of electron density variations along 
the different lines of sight in the interstellar medium. How- 
ever this is very unlikely, because the variations in DM caused 
by irregularities in the Galactic electron column density are at 
the level 5DM ^ 0.05 cm~^ pc over the small angle subtended 
by the pulsars (Nordgren, Cordes, & Terzian 1992), while the 
variations we observe are 10 times larger than this. 

We conclude that the most likely cause of the observed DM 
variations is the presence of a free electron plasma within the 
cluster, and estimate its density in the next section. 



Fig. 2. — Measured DM plotted against (^/^)obs ^™ ^^ millisecond pulsars. 
Most of the pulsars with negative (P/P)obs have significantly higher DMs than 
those with positive (P/P)obs ■ The curves show the expected observed variation 
of (P/P) with DM along lines of sight passing through the cluster core (dotted) 
and passing 2 pc from the cluster core {dashed). These are calculated using the 
cluster mass distiibution described in Freire et al. (2001), assuming that all the 
pulsars have the same intrinsic ((P/P)\nt), and that there is a uniform electron 
density of 0.067 cm"' throughout the cluster. The vertical and horizontal lines 
represent the fitted values of ((P/P)i„t) and DMc as described in § pi 



3. PLASMA DENSITY 

We derive the plasma density in the central regions of 47 Tuc 
using two different methods. In the first, and simpler method, 
the density n^ of a uniform plasma is obtained from the ratio 
of dispersion in measured DMs to the one-dimensional disper- 
sion in offsets in the plane of the sky from the cluster center 
(Table |I]): 



ctDM 0.13±0.04cm-3pc 



(t(9±D) 0.75±0.14pc 



: 0.17 ± 0.05 cm- 



(1) 



where we have used D = 5 kpc (Reid 1998 ), and assume a spher- 
ically symmetric distribution of pulsars. 

In the second method we calculate the positions of the pulsars 
along the line of sight from their period derivatives and compare 
these with the DMs. The measured values of P/P are biased by 
the intrinsic period derivative of each pulsar: with no a priori 
knowledge of {P/P)iai for individual pulsars, we assume that 
all have the average value {{P/P)mt)- Then, from the inferred 
accelerations a/c = {P/P)obs- {iP/P)mt) and a King model for 
the gravitat ional potential of the cluster (|King 1966|; see F reire 
et al. 2001), we compute the distance R of each pulsar along 
the line of sight from the plane of the sky containing the cen- 
ter of the cluster Here we are assuming that a King model 
for the potential computed with known cluster parameters, to- 
gether with a relatively small contribution from {P/P)int, pro- 
vides a good description for the observed values of P/P, which 
Freire et al. (2001) have shown to be the case. Using a uniform 
electron density n,,, we then calculate the incremental DM with 
respect to the DM to the center of the cluster, DMc, required 
to obtain the observed DM for each pulsar. In summary, our 



Freire et al. 



model for the DMs contains the free parameters DMc, He, and 
{(P/P)[nt)' ^nd is represented by 



VMi = DM, + n,Rii[P/PUsi- {[P/PU)), 



(2) 



the cloud of ionized gas at the center of 47 Tuc is fairly homo- 
geneous on r^ pc scales, and that the King model ( one simple 
choice amid many possible ones; see ^hinney 1993 ) provides a 
consistent description for the potential of the cluster 



where we have noted the explicit dependence of each pulsar's 
line-of-sight distance Rj on the inferred acceleration. There are 
often two distances at which a given a/c can be obtained (see 
Fig. g). We resolve this ambiguity by choosing the distance 
that produces a model dispersion measure VAi closest to the 
observed value DM. 

We determine the model parameters by minimizing 



fl: 



J2i^Mi-VMif 



(3) 



/=i 




and in order to obtain reliable uncertainty estimates in the de- 
rived parameters we implement this with a Monte Carlo pro- 
cedure. We do so by generating 10,000 data sets where both 
the observed pulsar DMs and relevant cluster parameters are 
chosen in a random fashion consisten t with their measured un- 
certainties (distance: 5.0±0.4kpc; Reid 1998; central stel- 
lar line-of-sight v elocity dispersion: 11.6± 1.4kms~'; M ey- 
lan & Mayor 1986; angular co re radius: 23."1 ± 1"4; H owell, 
Guhathakurta, & Gilliland 2000). 

Minimization of /i (eq. [ph) in the Monte Carlo procedure 
resulted in an average and rms for the model parameters^ of 
n, = 0.067 ± 0.015 cm-^ DM^ = 24.381 ± 0.009 cm'^pc and 
{(P/P\J = (0.48 ± 0.18) X 10"'^ s"'. The latter corresponds 
to an average "characteristic age" of {{P/2P)^^^) = 3 Gyr and an 
average magnetic field of « 3 x 10^ Gauss for these 15 mil- 
lisecond pulsars in 47 Tuc. The line-of-sight distances /?, and 
uncertainties resulting from the fitting procedure are displayed 
in Figure H and listed in Table n^. 

Each pulsar's {P/P)^^^^ is different from the average. In or- 
der to investigate the effect of this approximation in the eval- 
uation of our model parameters we repeated the Monte Carlo 
procedure described above, replacing {{P/P)^^^ with a value 
of {P/P)[ni for each pulsar drawn from a Gaussian distribu- 
tion of characteristic ages centered on 5 Gyr with an rms of 
2.5 Gyr, truncated at the low end by the limits derived by Freire 
et al. (2001). We repeated the procedure for a flat distribution 
of ages between 2 and 9 Gyr Both distributions were chosen 
to be approximately representative of the distribution of char- 
acteristic ages observed for ^ 30 millisecond pulsars in the 
Galactic disk. The results were similar in both cases, yield- 
ing parameters smaller by about 1 a compared to the original 
model represented by equation (0): n^ = 0.058;!;og[2 cm"^ and 
DMc = 24.37 ±0.02cm"^pc. Additionally, for three pulsars 
the inferred distances Ri or their uncertainties changed signif- 
icantly, owing to the prominence of the "second solution" in 
Ri{a) (cf. Fig. g), and are indicated in Table [T| 

Individual DMs are affected by fluctuations in the Galactic 
column density, as noted in § g. There is also some uncertainty 
regarding the detailed form of the potential used to derive Ri 
from the inferred accelerations. And each pulsar's character- 
istic age is different from the average value determined in the 
model. All these factors contribute to the scatter visible in Fig- 
ure g. Nevertheless, the small scatter about the linear relation 
between DM and R emphasizes that the observations are con- 
sistent with two of the central assumptions made in the fit: that 

^In general the parameters obtained in this procedure are not Gaussian distributed. Where the difference is small we approximate to average and rms values, and 
otherwise present the median and one-sided "34%" error bars. 



Fig. 3. — The measured DM plotted against the inferred line-of-sight dis- 
tance of 1 5 pulsars from the plane of the sky containing the cluster center A 
strong correlation is observed, for which the inferred free electron density is 
0.067 ± 0.015 cm-3. 



The values of n,, derived by the methods represented by equa- 
tions ([l|) and (g) agree at the 2 a level, although the value de- 
rived in equation ([l|) is nominally higher than the one derived 
in the Monte Carlo procedure. This is partly due to the appar- 
ently smaller effective radius of the pulsar distribution in the 
plane of the sky (1.9pc; Fig. n^), than along the line of sight 
(3.3 pc; Fig. ^. We attribute this apparent asphericity of the 
pulsar distribution to small-number statistics (confirmed with 
a KS test), and hereafter consider the pulsars to be distributed 
in a sphere with effective radius the geometrical mean of the 
two, rp = 2.5 pc. We also retain 0.067 ± 0.015 cm"-' as the most 
precise estimate for Ug. 

Eight pulsars with well-measured DM and {P/P)obs are also 
known in the globular cluster M15. The DM variations among 
these pulsars were interpreted by Anderson (1992) as arising 
from a gradient in the Galactic electron column density across 
the cluster. As an alternative interpretation, we find the data for 
the four millisecond pulsars (P < 10 ms) in M15 to be consis- 
tent with the existence of a homogeneous plasma with an elec- 
tron density of about 0.2 cm~^, slightly larger than in 47 Tuc. 
The other pulsars in M15 have longer periods and likely have 
large contaminating values of iP/P)inf 

4. DISCUSSION 

Using high precision pulsar timing applied to 15 pulsars ob- 
served in the globular cluster 47 Tuc, we have detected ion- 
ized intracluster gas with a density undetectable by other meth- 
ods. Our meas urement is consistent w ith previous upper lim- 
its for 47 Tuc (Hesser & Shawl 1977). The central electron 
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density derived, tie = 0.07 cm , compares to ^ 0.0007 cm for 
the ISM in the vicinity of 47 Tuc, at a z-height of 3.5 kpc, an d 
^ . 02 cm~^ in the solar neighborhood ( Taylor & Cordes 1993 ) . 
With a value of 0.07 cm""* for the mean electron density and 
assuming that each electron is accompanied by a single pro- 
ton, we obtain a plasma mass density p ~ 10"^^ g cm""*. Within 
the approximately spherical region of radius rp = 2.5 pc occu- 
pied by the pulsars this co rresponds to a total ma ss ^0.1 Mq. 
Recent ISO observations ( [Hop wood et al. 1999 ) put a strong 
upper limit of 2 x IQ~*Mq on the dust content in 47 Tuc, so 
that virtually any intracluster medium should be present as ion- 
ized or neutral gas. Moreover, we expect any gas present to 
be completely ionized: the detection of hot st ars in the cluster 
( |Pai-esce et al. 199l|;|Guhathakurtaet al. 1992|; O 'Connell et al. 
1997) indicates an intense background of UV photons, enough 
to ionize at least lO'^'^M© yr"'. This is more than enough to 
completely ionize the M ^ 10~^ Mc^ yr~' which i s lost by all the 



evolved stars in a typical cluster (Roberts 1986). We therefore 
conclude that the detected plasma represents the total intraclus- 
ter medium in the inner few pc. 

The inferred ^0.1 M0 is much less than the ^ 100 Mq of in- 
tracluster material expected to accu mulate within rp of the clus- 
ter core over a period of 10^-10^ yr ( Roberts 1986 ). What pro- 
cess is then responsible for ejecting most of the gas from this re- 
gion? Several different mechanisms have been proposed, most 
of wh ich invoke w inds suc h as those driven by ma in-sequence 
stars (ISmith 1999|), novae ( |Scott & Durisen 1978| ), M dwarfs 
(Coleman & Worden 1977), or pulsars. Indeed, Spergel (1991) 
proposed that 47 Tuc and other globular clusters should be de- 
void of any plasma, owing to the kinetic effects of strong winds 
from millisecond pulsars. Using the derived {{P/P^^^), the typ- 



ical spin-down luminosity (ex P/P^) of a millisecond pulsar in 
47 Tuc is w lO-'"^ erg s"' . Thus, the ^ 10^"* erg s"' needed to ex- 
pel 10"^ Mq yr"' from the clu ster's potential, which requires an 
escape velocity of 58 km s"' (Webbink 1985), can be provided 
by just ~ 0.5% of the spin-down l uminosity of the ^ 200 pul- 
sars believed to exist in the cluster (Camilo et al. 200C ). 



If we assume that mass loss within the central few pc of 
t he cluster is M ^ / x 10"^ Mq yr"', with / of order unity 



(Roberts 1986), then a steady state situation implies that the 
plasma is being expelled from the inner rp = 2.5 pc with ve- 
locity V = M/(47rr^p) ^ / x 80kms"', which is comparable to 
the escape velocity of the cluster. Also, the existence of this 
plasma, together with sensitive X-ray limits on accretion lumi- 
nosity from a central source, can be used to place an upper limit 
on the ma ss of a central black hole of about 100 Mq (G rindlay 
et al. 2001). 

To conclude, we have detected the long-sought intracluster 
medium in a globular cluster, 47 Tucanae. It is somewhat ironic 
that the very objects which allow us to detect the ionized gas 
might also be responsible for its low concentration in the clus- 
ter. 
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Table 1 
Parameters for 15 millisecond pulsars in 47 Tuc 



Pulsar 


(P/P)oh. 


0± 


DM 


R' 




(lO-'^s-i) 


(arcmin) 


(cm"-' pc) 


(pc) 


C 


-0.91 


1.21 


24.599(2) 


+2.6 ±0.6 


D 


-0.10 


0.68 


24.729(2) 


+5.2-]i 


E 


+2.74 


0.65 


24.230(2) 


-1.5±0.4 


F 


+2.42 


0.19 


24.379(5) 


_A 1 5+0.04 


G 


-1.09 


0.29 


24.441(5) 


,n 1 9+0.08 


H 


-0.09 


0.77 


24.36(3) 


+0.04 ±0.02 


I 


-1.36 


0.29 


24.42(2) 


+o.i5!»:i^ 


J 


-0.51 


1.00 


24.5867(7) 


+3.4!«:? 


L 


-2.85 


0.14 


24.38(5) 


+0.3!°:? 


M 


-1.08 


1.05 


24.42(2) 


+o.i2!«:«^ 


N 


-0.76 


0.49 


24.56(3) 


+7 Q^"^ 


O 


+1.10 


0.06 


24.368(8) 


-0.04 ±0.02 


Q 


+0.80 


0.98 


24.29(2) 


-0 n?+""i 


T 


+3.84 


0.34 


24.39(6) 


-0.4 ±0.2 


U 


+2.15 


0.94 


24.335(4) 


-1 2+0-3 
^ -^-0.2 



Note. — Uncertainties in {P/P)ohs and d± are less than one in 
the least-significant digits. The values of (P/P)obs listed have been 
obtained from the actually observed values after a corr ection for 



centrifugal and Galactic acceleration (Freire et al. 2001 ). The un- 
certainties in the last quoted digits of the DMs are given in paren- 
thesis. 

"The line-of-sight distances from the center of cluster (R > 
for the distant half) are model-dependent. We list values obtained 
with the model using a uniform characteristic age, represented by 
equation (0). For the following three pulsars, the distance or its un- 
certainty is different when obtained with the model where (P/P)int 
is drawn from a distribution (see § Bb: i^c = +0.2!^ j, Rm = +0. l^^ j, 

Ru=-O.UTd- 



